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Abstract 

We show for a large class of discrete Harper-like and continuous magnetic Schrodinger 
operators that their band edges are Lipschitz continuous with respect to the intensity of the 
£^ . external constant magnetic field. We generalize a result obtained by J. Bellissard in 1994, 

^ ' and give examples in favor of a recent conjecture of G. Nenciu. 



^ ! 1 Introduction and the main results 

Harper-like operators. Let T C IR 2 be a (possibly irregular) lattice which has the property that 
Oh there exists an injective map F : T n- Z 2 such that ^(7) — 7I < 1/2. The Hilbert space is / 2 (r). 

^ ch The elements of the canonical basis in l 2 (T) are denoted by {£ x }xgr, where <5 x (y) = 1 if y = x 

and zero otherwise. In the discrete case, to any bounded self-adjoint operator H E B(l 2 (T)) it cor- 
responds a bounded and symmetric kernel _ff(x, x') = (H5 x r,5 x ) = if(x',x). We will extensively 
use the Schur-Holmgren upper bound for the norm of a self-adjoint operator: 
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ff||< sup^|ff(x,x')|. (1.1) 



Denote by (x — x )" = [1 + (x - x ) 2 ]^, a > 0. We define C a to be the set of bounded and 
self-adjoint operators H € B(l 2 (T)) which have the property that their kernels obey a weighted 
' Schur-Holmgren type estimate: 

||.ff||c.> := sup y(x-x') Q |ff(x,x')| <0 o. (1.2) 



We also define the space H a which contains bounded and self-adjoint operators H which obey: 



><' \\B:\\h* ■■= sup V<x- X ') 2a |ff( X ,x')| 2 <oo. (1.3) 

' x'er 1 ' 1 



,xer J 

The flux of a unit magnetic field orthogonal to the plane through a triangle generated by x, x' 
and the origin is given by: 

<ys(x,x') := -i (xi x' 2 - x 2 x[) = -(p(x',x). (1.4) 
Note the important additive identity: 

<y5(x, y) + ip(y, x') = ip(x, x') + tp(x - y, y - x'), (1.5) 

\<p(x-y,y -x!)\ < ^l x -y| |y- x 'l- 

Let K 6 C°. Let its kernel be K fx, x'). We are interested in a family of Harper-like operators 
{Kb}beR given by the kernels e lb ^( x ' x 'iT(x, x'). Clearly, {Kb}beR C C°. The usual Harper operator 
lives in l 2 (Z 2 ), and its generating kernel has the form K(x, x') = k(x — x') where fc(x) equals 1 if 
|x| = 1, and otherwise. 

In Lemma [2.11 we will show that T~L a C C° if a > 1. Now here is the first main result of our 
paper: 
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Theorem 1.1. Let a > 3 and K £ % a . Construct the corresponding family of Harper-like 
operators {Kf,}beR- Then we have: 

i. The resolvent set p(Kb) is stable; more precisely, if dist(z, a(K{, )) > e > then there exist 
S > and rj > smc/i i/iai dist(z, o~(K{,)) > rj whenever \b — bo\ < S. 

ii. Define E+(b) :— sup o~(Kb) and E_(b) := inf o-(Kb). Then E± are Lipschitz functions ofb. 

iii. Let a > 4. j4ssMme t/iat -Kb /ias a gap m t/ie spectrum of the form (e_(&o), e+(&o))> where 
e ±(bo) S cr(-Kb ) are i/ie gap edges. Then as long as the gap is not closing by varying b in a closed 
interval I containing bo, the operator Kb will have a gap (e_ (b), e+(6)) whose edges are Lipschitz 
functions ofb on L. 

Remark. Denoting by 5b = b — bo, then according to our notations we have that Kb = {Kb ) 5b . 
It means that it is enough to prove spectral stability and Lipschitz properties near bo = 0. 

We can complicate the setting by allowing the generating kernel to depend on b. 

Corollary 1.2. Assume that the generating kernel K(x, x';6) obeys all the spatial localization 
conditions of Theorem ] 1. 11 uniformly in b G M. Moreover, assume that it also satisfies an extra 
condition: 

sup V|#(x,x';&)-#(x,x';& )| <C|b-6 |, |&-&o|<l. (1-6) 
x ' er xer 

Consider the family {Kb}be& generated by e 26 ^*'* ^-K"(x, x'; b). Then Theorem \l.l\ holds true for 
K b . 

Continuous Schrodinger operators. Let us consider the operator in L 2 (M. 2 ) 

H (b) := (p - fea) 2 + V, p = -iV x , a(x) = (-x 2 /2, x t /2), b e R. (1.7) 

where we assume that the scalar potential V is smooth and bounded together with all its derivatives 
on R 2 . This very strong condition is definitely not necessary for the result given below, but it 
simplifies the presentation. For the same reason we formulate the result only near bo = 0. 

Theorem 1.3. Assume that the spectrum of H{Q) has a finite and isolated spectral bando~o, where 
(To = [s— (0), s_|_(0)] . Then if \b\ is small enough, o~o will evolve into a still isolated spectral island 
o~b C o~(H(b)). Denote by s_ (6) := inf <r& and s + (b) := supcb. Then these edges are Lipschitz at 
b = 0, i.e. there exists a constant C such that \s±{b) — s±(0)| < C \b\. 

Remark. We do not exclude the appearance of gaps inside a\,. Moreover, the formulation 
of this result is slightly different from the one we gave in the discrete case. Here we look at the 
edges of a finite part of the spectrum, and not at the edges of a gap. In the discrete case both 
formulations are equivalent. However, our proof does not work in the continuous case if o~o is 
infinite. 

1.1 Previous results and open problems 

Spectrum stability is a fundamental issue in perturbation theory. It is well known that if W is 
relatively bounded to Ho, then the spectrum of H\ = Ho + XW is at a Hausdorff distance of 
order |A| from the spectrum of Ho- But this is in general not true for perturbations which are not 
relatively bounded. And the magnetic perturbation coming from a constant field is not relatively 
bounded, neither in the discrete nor in the continuous case. 

With the notable exception of a recent paper by Nenciu [5S] , all previous results on the discrete 
case we are aware of deal with the situation in which F = Z 2 and the generating kernel obeys 
K(x, x') = fc(x — x'), where k is sufficiently fast decaying at infinity. Maybe the first proof of 
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spectral stability of Harper operators is due to Elliott [TT] . The result is refined in where it is 
shown that the gap boundaries are ^-Holder continuous in b. Later results by Avron, van Mouche 
and Simon [2J, Helffer and Sjostrand [TBIE], and Haagerup and R0rdam [T3] pushed the exponent 
up to i . In fact they prove more, they show that the Hausdorff distance between spectra behaves 
like | — foo I ^ - These results are optimal in the sense that the Holder constant is independent of the 
length of the eventual gaps, and it is known that these gaps can close down precisely like \b — bo\ 2 
near rational values of 60 [HI HSj ■ Note that Nenciu [35] proves a similar result for a much larger 
class of Harper-like operators. Many other spectral properties of Harper operators can be found 
in a paper by Herrmann and Janssen |18j . 

In the continuous case, the stability of gaps was first shown by Avron and Simon [I], and 
Nenciu [37]. Nenciu's result implicitly gives a ^-Holder continuity in b for the Hausdorff distance 
between spectra. Then in [4i the Holder exponent of gap edges was pushed up to §. 

The first proof of Lipschitz continuity of gap edges for Harper-like operators was given by 
Bcllissard [3] (later on Kotani [21] extended his method to more general regular lattices and 
dimensions larger than two). The configuration space is V = 1? and the generating kernel is of 
the form K (x, x') = £;(x — x'; b), where fc(x; b) decays polynomially in |x| and is allowed to depend 
smoothly on b. This extra-dependence is not central for our discussion, so we will consider that k 
is b independent. Bellissard's innovative idea uses in an essential way that the Harper operators 
generated by translation invariant and fast decaying kernels fc(x — x') can be written as linear 
combinations of magnetic translations: 

K b = HlWbh), [Wi(7M(x) = e^ (x ^V(x - 7), W 6 (7)W 6 (7') = e ih ^'^Wb{l + 7'). 
7ez 2 

Bellissard's crucial observation was that the C* algebra Ab +s generated by {Wb +s{l)}-yez 2 is 
isomorphic with a sub-algebra of Ab <8> As which is generated by {W / h ( 7 ) ® Ws{"f)} ie z 2 - Thus 
one can construct an operator Kb +s which is isospectral with Kb a +s- The new operator lives in 
the space l 2 (J?) eg L 2 (R), and Kb = Kb £§> Id. It turns out that it is more convenient to study the 
spectral edges of the new operator. The reason is that the singularity induced by the magnetic 
perturbation is hidden in the extra-dimension. But the proof breaks down in case of irregular 
lattices or if the generating kernel K (x, x') is not just a function of x — x'. 

Coming back to our proof, its crucial ingredient consists in expressing the magnetic phases with 
the help of the heat kernel of a continuous Schrodinger operator, see (|5.8[) - (|5.12|) . Moreover, the 
proof in the discrete case also works for continuous kernels living on M 2 and not just on lattices. 
This is what we use in the last step of the proof of Theorem 11.31 dealing with continuous magnetic 
Schrodinger operators. 

A limitation of our method consists in the fact that the phases </?(x, x') are generated by 
a constant magnetic field. A more general discrete problem was formulated by Nenciu in 28 
where he proposed to replace the explicit formulas in (jl.4l) and (11.51) with more general real and 
antisymmetric phases obeying </>(x, x') = </>(x, x') = — </>(x',x) and 

\4>{x, y) + (f>(y, x') + 0(x', x)| < area A(x, y, x') 

where A(x, y, x') is the triangle generated by the three points. These phases appear very naturally 
in the continuous case, see [7j|5J[]J2[22l[23l[Ml[25l|26], where it is shown that if a(x) is the transverse 
gauge generated by a globally bounded magnetic field |6(x)| < 1, then 4>(x, x') can be chosen to 
be the path integral of a(x) on the segment linking x' with x. This is the same as the magnetic 
flux of b through the triangle generated by x, x' and the origin. 

Using a completely different proof method, Nenciu shows among other things in [28] that the 
gap edges are Lipschitz up to a logarithmic factor, and he conjectures that they are actually 
Lipschitz. His method relies on the theory of almost convex functions, and the result provided 
by this technique is optimal in the sense that it cannot be improved in order to get rid of the 
logarithm. A new idea would be necessary in order to prove Nenciu's Lipschitz conjecture. 

Our current paper supports this conjecture because it provides examples of phases not coming 
from a constant magnetic field which still generate Lipschitz gap edges. Let us show this here. 
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Consider an irregular lattice T C R 2 which is a local deformation of 1? , that is there exists a bi- 
jectivemapF : T — > I? such that \F{'y)— 7| < \- Define the phases ^(x, x') := <^(F _1 (x), F _1 (x')) 
where ip is given by (| 1 .4[) . 

Choose any self-adjoint operator K £ B(l 2 (Z 2 )) given by a kernel A(x, x') sufficiently fast 
decaying outside the diagonal. The same operator can be seen in B(l 2 (T)) given by if (7, 7') := 
K(F(-y),F (7')). Thus the operator K b generated by K b (x,x') := e l ^ (x ' x,) A'(x, x') is unitary 
equivalent with an operator in B(l 2 (T)) with a kernel 

X 6 (7, 7 ') :=e ib ^^K( in '). 

In this case, we know from Theorem 11.11 that the edges of the spectral gaps of AT, and thus Kb 
will have a Lipschitz behavior. But the general case remains open. 

2 Proof of Theorem 11.11 

This section is dedicated to the proof of our first theorem. Parts of this proof will be later on 
adapted to the continuous case in Theorem 11.31 

2.1 Proof of (i) 

Let us start by showing the existence of natural embeddings of C a 's in H a, s given by the following 
short lemma: 

Lemma 2.1. Let H £ T-L a with a > 1. Then H £ C 3 with (3 < a — 1. In particular, if a > 3 
then the kernel (x — x') 2 |if(x, x')| obeys a Schur- Holmgren estimate and thus defines a bounded 
operator. 

Proof. Choose some small enough e > such that a > f3 + 1 + e. We write: 
(x-x') /3 |-ff(x,x')| < (x-xV 1_e (x-x') Q |i/(x,x')| 
and see that the Cauchy-Schwarz inequality gives 

\\H\y<C a , ? \\H\\ na . (2.1) 

□ 

Another technical estimate to be proved in the Appendix claims that if H has a kernel which 
is localized near the diagonal, then the resolvent's kernel will also have such a localization. 

Proposition 2.2. Let H G C a , with a > 0. Let z £ p(H). Then for every < a' < a we have 
(H — z)~ l £ H a , and there exists a constant C independent of z such that 

\\(H - z)-^ <C(1 + ( {dist( J m}a+2 + dist(z ; g(g)) ) • (2-2) 

Remark. This proposition is related to what specialists in von Neumann algebras would call dual 
action, see [29] [301 [31] . Stronger localization results have been earlier obtained by Jaffard [20 , 
later generalized by Grochenig and Leinert [13] . We choose for completeness to give an elementary 
proof in the Appendix; our proof also highlights the uniformity in z £ p(H). 

Now let us start the proof of (i). Constants only depending on e will be named C e even though 
they might have different values. 

Remember that it is enough to prove the stability result near bo = 0. Let K £ H a with a > 3. 
Lemma [2.11 gives us some f3 > 2 such that K £ C 3 . Proposition 12.21 savs that (K — z)^ 1 £ *HP 
with some 2 < j3' < /3, while Lemma |2~T1 insures that there exists 7 > 1 such that (K — £ C 7 . 
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Denote by G(x, x'; z) the kernel of (K — z) 1 . From (|2.2[) and (|2.1j) we obtain a constant G e 
such that: 

sup V(x-x')|G(x,x';z)| < C e if dist(^, a(K)) > e. (2.3) 
x ' er xer 

Define the operator S b (z) to be the one corresponding to the kernel e lhip ( x ' x ^G(x, x'; z). Using the 
Schur-Holmgren criterion we can write 

\\S b (z)\\ < C e , beR, dist(z,a(K))>e. 

Using (| 1 . 5|) we can write: 

(K b -z)S b (z)=:l+T b {z), (2.4) 

where T b (z) is given by the kernel 

e *Mx,x') ^ ( e 'Mx-y,x'- y ) _ i)if( X) y ) G(y, x'; z). (2.5) 
yer 

Note that 

| e ^(x- y ,x'-y) _ !| < |^ (x _ y x / _ y) | < |l | x _ y | | y _ x /|. (2 6) 

Then for any / 6 ^ 2 (r) with compact support we can write: 

|T fc (z)/|(x) < |6|^|x-y| |Jf(x,y)| |y-x'| |G(y,x';z)| |/(x')| (2.7) 
yer 

and after applying the Schur-Holmgren criterion we get: 

\\T b (z)\\<\b\\\K\\ e ,\\(K-z)- 1 \\ c ,<\b\C e . 

Thus if \b\ is small enough, ||Tb(z)|| < 1/2 whenever dist(z, a(K )) > e. Now if Imz ^ we know 
that Kb — z is invertible, and from (|2.4[) we conclude that there exists a constant C e such that 

{K b - z)- 1 = S b (z) (l+T b (z)y\ 
1 1 (if 6 - z) _1 || < G e whenever |6| < b e and dist(z,a(K)) > e, (2.8) 

uniformly in the imaginary part of z. This means that dist(z, a(K b )) > ^- > whenever |6| < b t 
and dist(z, CT(if)) > e, and the proof of (i) is over. □ 

2.2 Proof of (it) 

As before, we only need to consider bo — 0. We give the proof just for the upper spectral limit 
since the argument for _E„ is similar. 

2.2.1 Reduction to localized operators 

We start with an abstract lemma. 

Lemma 2.3. Let M(b) and N(b) be two families of bounded and self-adjoint operators on some 
Hilbert space U, such that \\M{b) - N{b)\\ < C \b\ if \b\ < 1. Then: 

\ sup a(M(b))- sup a(N(b))\<\\M(b)-N(b)\\<C\b\, \b\<l, (2.9) 

and a similar estimate holds for the infimum of their spectra. In particular, if sup a(N(b)) is 
Lipschitz at b = then the same is true for suptr(M(6)). 

Remark. Note the important thing that we do not require from M(b) and N(b) to converge in 
norm to M(0) = N(0) when b tends to zero. 
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Proof. For every if) G "H with H^H = 1 we can write 

(M(b)i>, < (N(b)iP,i>) + \\M(b) - N(b)\\ < sup a(N(b)) + \\M(b) - N(B)\\ 

which means that sup a{M{b))- sup cr(N(b)) < \ \M(b)~ N(b)\\. By interchanging M(b) with N(b) 
we obtain the inequality: 

| sup a{M(b)) - supcr(7V(&))| < \\M(b) - iV(6)||. (2.10) 

A similar argument shows the same estimate for the infimum of the spectra. Regarding the Lips- 
chitz property, we use that supcr(M(0)) = sup cr(./V(0)) and then we apply the triangle inequality: 

|supa(M(6)) -sup<r(M(0))| < | sup<r(iV(6)) - supcr(iV(0))| + \\M(b) ~ N(b)\\ < C \b\. (2.11) 

□ 



Getting back to our theorem, we now want to reduce the problem to operators with kernels 
supported near the diagonal. Denote by x the characteristic function of the interval [0, 1]. Denote 
by Kb the operator given by the kernel Kb(x.,x.') := ^(Vfolx — x'|)X(x,x') and by Kb the operator 
given by K b (x,x') := x(v / 5|x - x'De^^i^x, x'). 

Since K £ H a with a > 3, according to Lemma [2~T1 we have the bound: 

(2.12) 



sup^(x-x') 2 |X(x,x')| = ||^|b < 

v 'czr — 



x'er 



xer 



Via the Schur-Holmgren criterion we obtain: 

maxiWK - K b \\, \\K b -kb\\} < sup V fl-x(V6|x-x'|)l |#(x,x')| < \b\ \\K\\ C * 



(2.13) 



xer 



Using Lemma 12.31 for the pair K and Kb we obtain |£ + (0) — sup(cr(i^£,)) | < \b\ \\K\\ C 2. The 
same lemma for the pair Kb and Kb gives \E + (b) — sup(er(i^b))| < \b\ ||X|| C 2. Then the triangle 
inequality leads to: 

\E+(b) - £+(0)| < 2\b\ \\K\\ C 2 + | sn V (a(Kb)) - sup(a(K b ))\. (2.14) 
Thus we have reduced the problem to the study of the spectral edges of Kb and K b . 
2.2.2 Study of the operators with cut-off 

Clearly, K b (x.,x.') = e ib ^ x ' x '^K b {x., x'). Without loss, assume that b > 0. Take if) £ l 2 (T) with 
compact support and compute (use (|5 . 10[) in the second equality): 

(K b i,,tP)= Y, e ibv(x ' x,) ^(x,x')^(x')?R 



x,x'er 



x,x'er 

Now denote by A b (t) the operator with kernel 

6|x-x'i 2 "1 



6|x — x 



'12 



4tanh(26t) 



G 6 (x,y;t)G&(y,x';t). (2.15) 



A b (x, x'; t) := K b (x, x') exp 



4tanh(2W) 



if (x, x') exp 



&|x - X' 



'12 



4tanh(26i) 



;(Vb\x-x'\ 
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The crucial observation is that equation (|2.15[) leads to: 



(K b iP, i>) = / dy(A b (t)G b (y, •; t)i/>, G b (y, •; t)^> 



4tt sinh(26i) 



R 2 



< sup a(A b (t)) 



sup a{A b (t)) 



4tt sinh(26i) 
b 

4tt sinh(26i) 



R 2 



dy||G 6 (y,-;^ll : 



R 2 



dyJ2\Gb(y,^t)\ 2 \^)f 



xer 



= supa(A(*)) 1 1 VI I 2 



(2.16) 



where in the last line we used (|5.12[) . It means that supa(K b ) < sup(a(A b (t))) for all t. Now let 
us show that the operator A b (t) — K b has a norm proportional with b if t is large enough (say 
t = 6 _1 ). Indeed, we can write 



lA^x';^ 1 )-^*,*')! < |if(x,x')|x(V5|x-x'|) (exp 
<|^(x,x')|x(V5|x-x'|) 



6|x — x 



'12 



4tanh(2) 



1 



„,6|x-x'| 2 r&|x-x'| 2 



exp 



(2.17) 



(2.18) 



4tanh(2) [4tanh(2) 
and on the support of \ we can bound the above difference with: 

\A b (x, x'; b- 1 ) - K h (x, x')| < const b |x - x'| 2 | J ftT(x, x') 
The right hand side defines an operator whose norm behaves like b. Thus (|2~T6|) and ([2~T8l) imply: 

sup a(K b ) < sup aiAbib- 1 )) and WA^ 1 ) - K b \\ < C b. (2.19) 
Using (|2.10p for the pair Ab(6 _1 ) and K b we arrive at: 

sup a(K b ) < sup cr(Jf 6 ) + C 6. (2.20) 

We now want to change places between K b and K b in the above inequality, which would lead 
to sup a(K b ) < sup a(K b ) + C b and thus: 

| sup a(K b ) — sup a{K b )\ < C b, 

which together with (|2.14[) would imply: 

\E+(b) - E+{0)\ < Cb, b>0. 

The key step in the proof of (g^P was (|2TT5|) . Since K b (x,x') = e- ib ^ x < x '^ K b (x, x') we can 
write (use (|5. in the second line): 



(K 6 V,^)= J] e -^( x ' x '^ b (x,x')V(x')^(x) 



x,x'er 



/v-^ , — ^^4-7r sinh(26i) ~ , ,. 
dy V>(x'M*) r^^(x,x')exp 

- 12 x,x'er 



6|x — x 



'12 



4tanh(26i) 



G 6 (x , ,y;t)G 6 (y,x;t). (2.21) 



Now everything will work as before, because the phase e 16 ^*'* ) changes neither the localization 
nor the C 2 norm of the operators. The proof for the upper spectral edges is over. 

The proof for the lower spectral edges is based on an estimate which is very similar with (|2.16p . 
in which we reverse the inequality and show that inf <r(K b ) > ini a(A b (t)) for all t. We give no 
further details. 
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2.3 Proof of (m) 

The idea is to reduce the problem to the previous case. Again it is enough to consider bo = 
and b > small enough. Assume that K has a gap in its spectrum of the form (e_,e+), with 
e± £ o~(K). Then due to (i) we know that if b is small enough the gap will survive: we can choose 
a positively oriented circle L in the complex plane containing := o~(K b ) n [e + (b),oo) such 

that 

dist(z, a (Kb)) > rj > whenever z £ L and < b < b„. 

The orthogonal projector P b corresponding to £+(&) can be written as a Riesz integral and we 
have: 



J(K b -z)- l dz, K b P b = ^-Jz(K b -z)- 1 dz, b>0. (2.22) 



If we consider KbPb as an operator living on the whole space l 2 (T), then its spectrum is given 
by the union {0} U E+(6). If we choose A := 1 + sup cr(K), then for b small enough the operator 
D b := K b P b - XPb will have inf a(D b ) = e+(b) - A < -1/2. Thus e+(b) = A + inf a(D b ), hence 
e+(6) is Lipschitz at b = if inf a(Db) has the same property. This is what we prove next: 

Lemma 2.4. Let Db = KbPb — XP b wzt/i A := 1 + sup o~(K). Then there exists b\ > small enough 
and a constant C > smc/i i/iaf /or every < 6 < &i we Ziawe | inf a(D b ) — inf <r(Do)| < C o. 

Proof. Remember that we imposed a > 4. We have that ||-?fb||-H<* = ||if||^a < oo for all 6. 
According to Lemma |2~T1 there exists j3 > 3 such that HifjUcP = ||if||c^ < oo. Then if b is smaller 
than some constant only depending on L, Proposition 12.21 tells us that (Kb — z)^ 1 £ Hr for some 
3 < P' < f3, for all z £ L and sup 2ei \ \(K b - z) _1 || H ^' < C. Thus both P b and D b belong to V. 13 ' 
with (3' > 3 if b is small enough. More precisely, there exists 62 > sufficiently small such that 

max{\\P b \\ nfi ,,\\D b \\ n e,}<C, < b < b 2 . (2.23) 

If G(x, x'; z) is the integral kernel of (K — z)^ 1 , then we introduced at point (i) the operator S b (z) 
given by the kernel e lb¥> ^ x,x -*G(x, x'; z). Using (|2.8I) we can write: 

sup 11(^6 - z)- 1 - S b (z)\\ <Cb, (2.24) 
provided 6 is small enough. Denoting by Dq the operator given by the integral kernel 

A)(x,x') := ^^(z-A)G(x,x';z)dz 

and by (Do) b the operator generated by e tbLp ^ x ' x 'Dq(x, x'), then using (|2.24[) we arrive at the 
estimate: 

H-D& - (£>o)b|| < G 6 whenever < b < b 2 . (2.25) 

It follows from Lemma l2~3"1 that inf o-(D b ) is Lipschitz at 6 = if inf a((Do)b) has the same property. 
But for the operator (-Do)fe we can apply point (ii), and the proof is over. □ 

3 Proof of Corollary [Q 

In order to keep the notation simple, we will only consider bo = 0. Here the generating kernel 
K(x, x'; 6) depends on 6 and (|1.6[) at 6 = reads as: 

sup V|i^(x,x';o)- A^(x,x';0)| < G |6|, |6| < 1. (3.1) 
x ' er xer 
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Let us introduce the family Kb where their kernels are given by e ib ^ x ' x ')A(x,x';0). Clearly, 
Kq = A'o- Moreover, (|3.1j) implies that \\Kb — K b \\ < C \b\ around b — 0. We know that Theorem 
11,11 (it) applies for Kb around b = 0, so the only thing we have left is to extend it to Kb- From 
Lemma 12.31 we immediately conclude that sup <y{Kb) an d inf cr(Ab) are Lipschitz at b = 0. 

The spectral stability of Kb can be shown with the same strategy as the one one used in (|2.3p ~ 
(I2.8[) . The operator Sb(z) must be constructed starting from the kernel of {Kq — z) _1 which gets 
multiplied with the phase e lbv( - x ' x When we act with Kb — z on Sb(z) as in (|2 .4[) . we obtain an 
extra term which enters in T b (z), which is (Kb — Kb)Sb(z). This error is again proportional with 
|6| if z is at some distance from the spectrum of Kq. Thus (|2.8[) holds again. 

For the case of gaps, the proof is identical with the case independent of b. □ 

4 Proof of Theorem 11.3 

There are important similarities between the proof strategies in the discrete and continuous cases. 
Although the stability of the resolvent set of H(b) is known, we will sketch a short proof which 
will also provide some ingredients for the proof of the Lipschitz behavior of the band edges. 

4.1 Stability of gaps 

Assume that M C p(H(0)) is a compact set and dist(M, a(H(0))) > 0. The resolvent (H(O)-zy 1 
is an integral operator given by an integral kernel Qo(x, x'; z). 

The singularities of Qq(x, x'; z) are the same as in the case of the free Laplacean and there 
exists some S > and Cm < oo such that uniformly in x / x' [321 HQ] : 

sup |Q (x,x';z)| < Cm (1 + | ln(|x - x'|)|) e - A >- x 'l, 

zeM 

sup |V x Q (x,x';z)| < C M U + r^—n) e - ' 1 *-*' 1 . (4.1) 
zeM \ \x-x'\J 

In particular we have the following Schur-Holmgren type property: 

sup sup / \Q (x, x'; z)\dx < C(M) < oo. (4.2) 

This allows us to define for every z € M a bounded operator Sb(z) whose integral kernel is given 
by: 

S b (x,x';z) :=e' ib ^ x < x )Qo(x,x';z), sup ||S 6 (z)|| < C{M) < oo. (4.3) 

zeM 

Define Tb(z) to be the operator with the integral kernel: 

T b (x, x'; z) := be lbv ^{2ta(x - x')V x Q (x, x'; z) + 6|a(x - x')| 2 g (x, x'; *)}. (4.4) 

The kernel Tb(x, x'; z) is bounded because the inequality |a(x — x')| < |x — x'| compensates the 
local singularities of V x <5o( x , x'; z) and <5o( x , x'; z) when |x — x'| is small, while when |x — x'| is 
large we have the exponential decay which comes into play. In fact, using (I4.ip we see that the 
kernel T b (x, x'\ z) obeys a Schur-Holmgren estimate. We get: 

sup ||T 6 (z)|| < C{M) \b\, \b\<l. (4.5) 

zeM 

Note the important identity valid on Schwartz functions: 

{-iV x - ba(x)}e tbv{x > x "> = e^ (x < x,) {-iV x - ba(x - x')}. (4.6) 
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Let us note that S b {z) leaves the Schwartz space invariant and for such two functions / and g we 
have (using (|4.6|1 ): 

({(p-b a ) 2 + V-z}S b (z)f,g) (4.7) 

= 11 e lMx ' x ' ) ({-«V x -6a(x-x')} 2 + V(x)-z)Q (x,x';z)/(x')^W ( ixdx' 
Jr 2 Jr 2 

= (f,g) + (T b (z)f,g). 

The operator H (b) is essentially self-adjoint on the Schwartz space, and after a density argument we 
conclude that the range of Sb(z) is contained in the domain of H (b) and (H(b)—z)Sb(z) = 1+T b (z). 
Now there exists b\ > small enough such that if |6| < b\ we have sup zgM ||Tf,(z)|| < 1/2 (see 
(|4.5|) V Then after a standard argument we conclude 

(tf(fc)-z)- 1 =S b (z)(l + T b (z))-\ sup \\(H(b) - z)-^ <C M , \b\<h. (4.8) 

zeM 

This means that the gaps in the spectrum of H (0) are preserved. In particular, for every e > 
there exists 62(e) > such that: 

s_(Q) - e < s_(6) < s + (b) < s+(0) + e whenever |6| < b 2 (e). (4.9) 

Choose a positively oriented circle L isolated from a(H(0)) such that L completely contains the 
finite band Uq. Then if \b\ is small enough L will completely contain a b and remain separated from 
a(H(b)). 

4.2 The reduction to Harper-like operators 

As in the discrete case, we construct the Riesz integrals 

Pb -ir- [ (H(b) - z)~ x dz, K(b) := H(b)P b = / z{H(b) - z^dz. 

The operator H(b)P b seen in the whole space L 2 (R 2 ) will have the spectrum a b U {0}. Fix 
A+ := 1 - s_(0). If \b\ < 6 2 (i) (see gl|)), then we know that s + (b) + X+ > s_(6) + A+ > \ > 0. 
It means that 

*+(&) + A+ = sup a{H(b)P b + X+P b }. 
Similarly, choosing A_ := —1 — s+(0) we have s_(6) + A_ < — ^ < hence 

s_(6) + A_ = inf a{H{b)P b + A_P fc }. 

In other words, the band edges s± will be Lipschitz at b = if the spectral edges of the operators 

K ± (b) := H{b)P b + X ± P b = / (z + \±)(H(b) - z)~ x dz 

have the same property. Note that the operator if±(0) has an integral kernel given by: 

X ± (0)(x, X ') = ±-J^ z + A±)Q (x, x'; «)dz, |^±(0)(x, x')| < Ce- 5 I X - X 'l , (4.10) 

where the local singularity at x = x' dissapears due to the integral with respect to z. 
Now using (|4~8)l . (|43)l and (|4~5j) we have: 



K -M-^~ I A I.S',,(:V/: 



<C|6|. (4.11) 
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According to Lemma 12. 31 the spectral edges of K± (b) are Lipschitz at b — if the same property 
is true for 

(K±(0)) b :=±J (z + \ ± )S b (z)dz. 

This notation wants to highlight the fact that (K±(0)) b is given by the integral kernel 

(K ± (0)) b (x, x') := e^( x ' x ')/f ± (0)(x, x'). 

At this point we are in a situation which is completely similar to the discrete case, with the 
difference that the Hilbert space is L 2 (R 2 ) and the sums over T have to be replaced by integrals. 
The unperturbed kernel K±(0)(x.,x') has an exponential localization. 

We can mimic the proof of Theorem ll.il (ii) and conclude that the spectral edges of (K±(0)) b 
are Lipschitz at b = 0, and we are done. □ 

5 Appendix 

5.1 Proof of Proposition 12.21 

Denote by G(x, x'; z) the integral kernel of (H — z)^ 1 . If a 1 = we have 

J2 |G(x,x';z)| 2 = \\(H- z^MI 2 < jrr- ( \ m U2 

^ {dlSt(z,cr(ff))} 2 

uniformly in x', an estimate which is in fact much better than (|2.2[) . So from now on we may 
assume that < a' < a. 

For k £ M 2 define the unitary multiplication operator U k by (U k f)(x) = e x /(x). Define 
the family of isospectral operators H k = U k HU£, with integral kernels given by Hk(x,x') = 
e jk-(x-x )/f(x ;X '). We need the following technical result: 

Lemma 5.1. Let H be an element ofC a . Let n be the integer part of a. Then the mapping 

R 2 9 ff k £ B(l 2 (T)) 

is n times continuously differentiable in the norm topology. Moreover, any n 'th order mixed partial 
derivative of is a — n Holder continuous at k = in the norm topology. 

Proof. Assume that k = (fci,^). The integral kernel of is e 4k ' x_x '-ff (x, x'). Let n be 
the integer part of a. Then is n times differentiable in the norm topology with respect to 
kj, j £ {1,2}, and its n'th mixed partial derivative d™ i d^~ m Hy i is given by the integral kernel 
i n (xi — x' 1 ) m (x2 — x' 2 ) n ^ m e lk ■( x_x )iJ(x, x'). This integral kernel defines a bounded operator 
because \(xi — x' 1 ) m (x2 — x 2 ) n ~ m \ < (x — x')" and then we can use fj 1 . 2|) . 

For the Holder continuity statement, we use the estimate |e lk '( x ~ x > — 1| < 2 1 ~ 13 |kp|x — x'|^ 
which holds for every < (3 < 1. 

□ 

Now let z £ p(H). Denote by Gk(x, x'; z) the integral kernel of (7Jk~ Due to the identity 
U k (H - z)- 1 ^ = (H k - z)- 1 we have: 

G k (x, x'; z) = e 4k - (x ~ x,) G(x, x'; z). (5.1) 

Let us denote by n the integer part of a. We can suppose that n > 1 since the case < a < 1 is 
covered by the argument below. 
From the identity 

{H V - z)- 1 - (H k - z)- 1 = -{H V - z)- 1 ^, - H k ](H k - z)- 1 (5.2) 
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and from Lemma |5. II we conclude that the map 

R 2 3k^ (ffk-z)- 1 G B(l 2 (T)) 
is continuous in the norm topology, and also differentiable. We have: 

D k (H k - zy 1 = -(H k - z)- x [D k H k ](H k - z)~ x . (5.3) 
Using this identity at k = in (|5 . 1 1) leads to: 

(x - x')G(x, x'; z) = -((H - zr^AcHkWff - z)' 1 ^, ix> 

which gives: 

1 1 



ii(tf-zni H1 <c(i+iitfii cl ) 



dist(z,cr(£f)) 2 dist(z, cr(H)) 



This is true because we have the pointwise bound 

(x - x') 2Q < (1 + la* - 4| + \x 2 - 4|) 2q < (3max{l, |x x - x'J, |x 2 - x' 2 |) 2q 

2 

<3 2a +^3 2a |^-x;| 2a . (5.4) 
i=i 

By induction we obtain the following rough estimate: 

ma - < c (i + ii^iisj ( dist( J (H)) „ t , + sj^jgjj) • m 

Now let us assume that n < a < n + 1. The integral kernel of the n'th partial derivative of 
(ifk — z) _1 with respect to fci is given by i™ e lk ( x - x _ x[) n G(x, x'; z). Moreover, using (15.31) 
and Lemma [5Tl we conclude that the operator d£ (i?k — is a — n Holder continuous at k = 0. 
Let k = (fci,0). We also have the identity: 

in(e iMxi-«i) _ _ x'^G^x'; z) = (K(ff fel - z)- 1 - dl(H kl - z)-\ =0 ]S x ,,6 x ). 

If | fci| < 1 the following norm estimate holds true according to Lemma \5. II 

UK (H„, - z)-' - dKH,,, -^-'Ik-olH 

< cw~(i + IfflK 1 ) ( dist( J (g)) „« + di5t( ^ (fl)) ) ■ m 

Choose n < a' < a < n + Then the following integral converges in norm and defines a bounded 
operator: 

■/O Ki 

Its integral kernel is given by 

G(x,x';z):- l "( 2 ; 1 -x 1 )"G(x,x';z) / ( e *(— ^ - 1)^. 

Jo »i 

Assuming without loss of generality that x\ — x\ ^ 0, and by a change of variable s = k\ |xi — x[\ 
we obtain: 

G(x,x';z) = |xx -xi| Q _n (*i -4rG(x )X ';^) / - TF _ r _i»( e » i) - l)da. 
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Notice that the above integral only has two possible values C± both different from zero, depending 
on the sign of x\ — x\. Since G(x, x'; z) = (H5 X >,5 X ) = C±(x\, x[) \x\— x[ \ a '~ n {xi~ x' 1 ) n G(x, x'; z) 
with \C±(x\, x[)\ > G it follows that 

sup £ |a: a - x[\ 2a ' |G(x,x'; z)\ 2 < C- 2 \\H\\ 2 . 
x ' er xer 

This argument can be repeated for the other coordinate and bound the I 2 norm of ( — x') Q G(-, x'; z) 
using (15.41) . The proof of Proposition [2^21 is over. 



5.2 A few identities from the continuous case 



We list here a few well known facts about the continuous two dimensional magnetic Schrodinger 
operator with constant magnetic field equal to b in L 2 (M. 2 ): 



H b = (p-fca(x)) 2 , p = -zV x , a(x) = (-x 2 /2,a; 1 /2). 



(5.7) 



The integral kernel of the semi-group e tHb is denoted with Gfc(x, x';£) and is given by the 
following explicit formula: 



G b (x,x;t) = e 



'• f\ = P ib v(x,x') 



■ exp 



6|x-: 



4tanh(6i) 



47r smh(bt) 

The semigroup property insures the following identity: 

G b (x, x'; 2t) = [ G b (x, y; t)G b (y, x';t)dy 

iR 2 



=: e ib¥, ( x ' x ')Gf,(x,x';t). 



Then we can write: 

ibtp(x,x'[ 



— f G b (x,y;i)G,(y,x';i)dy 

x, x';2£) Jr2 

fe|x-x'| 2 



G 6 (x,x';2t) 
4tt sinh(2itf) 



■ exp 



4tanh(2W) 



G b (x,y;t)G b (y,x.';t)dy. 



Taking the complex conjugation in both sides gives: 



-ibip(x.,x.') 



47rsinh(26i) 



■ exp 



&|x - x 



'12 



4tanh(26t) 



R 2 



Again the semi-group property gives that: 
b 



47rsinh(2ta) 
which is clearly x independent. 



G b (x, x; 2t) = f G b (x, y; t)G b (y, x; t)dy = f \G b (y, x; t)\ 2 dy 



(5.8) 



(5.9) 



(5.10) 



G fe (y,x;t)G b (x',y;i)dy. (5.11) 



(5.12) 
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